Thermal Casimir effect between random layered dielectrics 
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We study the thermal Casimir effect between two thick slabs composed of plane-parallel layers 
of random dielectric materials interacting across an intervening homogeneous dielectric. It is found 
that the effective interaction at long distances is self averaging and is given by a description in 
terms of effective dielectric functions. The behavior at short distances becomes random (sample 
dependent) and is dominated by the local values of the dielectric function proximal to each other 
across the dielectrically homogeneous slab. 

PACS numbers: 05.40.-a, 77.22.-d 



Systems with spatially varying dielectric functions ex- 
hibit effective van der Waals interactions arising from 
the interaction between fluctuating dipolcs in the system 
[TJ [2]. These fluctuation interactions have two distinct 
components: (i) a classical or thermal component due to 
the zero frequency response of the dipoles and (ii) a quan- 
tum component due to the non zero frequency/quantum 
response of the dipoles. Despite the clear physical dif- 
ferences in these contributions, the mathematical com- 
putation of the corresponding interaction is almost iden- 
tical and boils down to the computation of an appro- 
priate functional determinant. The full theory taking 
into account both of these component interactions is 
the celebrated Lifshitz theory of van der Waals inter- 
actions [3], based on boundary conditions imposed on 
the electromagnetic field at the bounding surfaces and 
the fluctuation-dissipation theorem for the electromag- 
netic potential operators. The original Casimir interac- 
tion [4] is obtained in the limit of zero temperature and 
ideally polarizable bounding surfaces. At non-zero tem- 
perature the contribution of the zero frequency modes to 
the Lifshitz theory yields the classical thermal Casimir 
effect which is due to the non-retarded van der Waals 
interactions. The major mathematical problems in the 
computation of Casimir type interactions (setting aside 
the experimental and theoretical challenges to determine 
the correct dielectric behavior) are (i) the application of 
the Lifshitz approach to non-trivial geometries and (ii) 
taking into account local inhomogeneities in the dielec- 
tric properties of the media, always present in realistic 
systems. In this paper we will address the latter. 

We consider the thermal Casimir interaction for the 
case where the local dielectric function is a random vari- 
able in the transverse direction. Specifically we will con- 
sider the interaction between two thick parallel dielectric 



slabs, separated by a homogenous dielectric medium, see 
Fig. ([I]). The thickness of both disordered dielectric 
slabs are L\ and L2 respectively and their separation is 
denoted by I. In what follows we will study the limit of 
infinite slabs i.e. L\, L 2 — > 00. The dielectric response 
within the two slabs is constant in the planes perpendic- 
ular to the slab normal, but varies in the direction of the 
surface normal. It is well known that this problem can 
be solved in the case where the dielectric constants of the 
slabs do not vary [5] and the result can be tentatively ap- 
plied to the case of fluctuating dielectric functions via an 
effective medium theory which consists of replacing the 
fluctuating dielectric functions by an effective (spatially 
constant within each of the slabs) dielectric tensor. The 
most commonly used approximation is that where the lo- 
cal dielectric tensor is replaced by the effective dielectric 
tensor [3H], i-e. 

eii(x)^eg>, (1) 

with the bulk dielectric tensor defined via e\j{Ej) = 
(eijEj). The use of the effective dielectric function is 
not easily justifiable mathematically as an approxima- 
tion, although physically the effective dielectric function 
clearly does capture the bulk response to constant electric 
fields. We shall see that, for the random layered dielectric 
model studied here, the effective dielectric constant ap- 
proximation of Eq. (JT|) does in fact give the correct value 
of the thermal Casimir interaction when the two slabs 
are widely separated. This can be expected on phys- 
ical grounds since the fluctuating electromagnetic field 
modes with smallest wave- vectors (corresponding to vari- 
ations on large scales) dominate the Casimir interaction 
for large inter-slab separation. The dielectric response of 
the material to a constant electric field is given by the 
effective dielectric constant and if the wave- vector depen- 
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FIG. 1: A schematic presentation of the model. Two finite 
slabs, (1) and (2), with disordered plane-parallel dielectric 
layers interacting across a dielectrically homogeneous slab of 
thickness £. z axis is perpendicular to the plane of the slabs. 



dent response is suitably analytic near k — we expect 
that $Ck)~eg>(Q)=$ for |k|«l. 

In this letter we introduce a path integral formalism 
to compute the thermal Casimir free energy between two 
semi-infinite dielectric slabs which are composed of layers 
with varying dielectric function. Our formulation allows 
us to show rigorously that for large inter-slab separations 
the leading order contribution to the interaction is self 
averaging and is equivalent to that obtained by replac- 
ing each slab with a homogeneous (though non-isotropic 
medium) with a dielectric tensor equal to the effective 
(bulk) dielectric tensor of the disordered medium. The 
short distance behavior of the interaction is random and 
we show, as would be expected on physical grounds, that 
it is dominated by the precise value of the dielectric con- 
stants at the two opposing slab faces. 

The Hamiltonian for the zero frequency fluctuations of 
the electrostatic field in a dielectric medium is given by 
the classical electromagnetic field energy 



H[cf>] = i / dx e(x) (V0(x)r 



(2) 



and the corresponding partition function is given by the 



functional integral Z = J d[4>] exp(— j3H[4>]). Differences 
in dielectric functions lead to the thermal Casimir effect. 
Here we will consider layered systems where the dielectric 
function e depends only on the z direction e(x) = e(z). 
If we express the field <p in terms of its Fourier modes 
in the plane perpendicular to z, and we take the area 
perpendicular to z as A, with wave- vector k = (k x ,k v ), 
then the Hamiltonian can be written as H = ^ k with 



dze(z) 



#(*,k) 



dz 



k 2 |0(z,k) 



(3) 



A direct consequence of this decomposition is that the 
partition function can be expressed as a sum over the 
partition functions of individual modes Z^ as ln(Z) = 
J2k m (^k) where 



Z k = J d[X]cxp \-~J 



dz e(z) 



dXY 

~dz~) 



k 2 X 2 



(4) 

Here k = |k| and we have taken into account that the 
field is real. 

The problem of computing the interaction between 
slabs composed of layers of finite thickness can be stud- 
ied using a transfer matrix method [5] . However we will 
use a method based on the Feynman path integral which 
is particularly well suited to the study of systems where 
the dielectric function can vary continuously [6j. If we 
specify the starting and ending points of the above path 
integral, we see that it has to be of a harmonic oscillator 
form defined by 



dX 
dz 



■co 2 X 2 



K{x,y\z) = j d[X]cxp ^-^ J dzM(z) 

(5) 

which can be computed using the generalized Pauli - van 
Vleck formula [6j [7] telling us that K must have the gen- 
eral form 



K(x,y;z) 



exp 



(6) 

We may now write down evolution equations for the coef- 
ficients a%, a/ and b using the Markovian property of the 
path integral, K(x, y; z+z') = J dw K(x, w; z)K(w, y; z') 
[51 [7], which can also be used to prove the generalized 
Pauli - van Vleck formula. We obtain the evolution equa- 
tions 



ddi 

dz 

da f 
dz 



db 

dz 



ba f 
M ' 



mJ 2 - -J- 

M 



(7) 
(8) 



We thus find that the ^-dependent part of the free energy 
of the mode k (up to a bulk term which can be subtracted 
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off to get the interaction energy) is given by 



{af{k)-e k){a>;>{k)-e a k) 



(2), 



In 1 



-2kt 



{af{k) + e Q k){a^{k)+e k) 



(2) 



and the total I dependent free energy is F 



(9) 

(12') l l 

Here ' (k) are the solutions to Eq. (8 1 evaluated at 
the opposing faces of each slab (1) and (ZJ respectively. 

In order to evaluate the integrals of a,j ' (k), one first 
has to solve equations of motion Eqs. ^ to get the z 
dependence of O/ffe, z) and then proceed to the integrals 
that enter Eq. The evolution equation for a/(fc) for 
cither slab can be read off from Eq. (|8|) and is given by 



dcif(k, z) 
dz 



(10) 



An appropriate Hopf-Cole transformation [7J shows 
this formalism to be equivalent to the transfer matrix 
method [5] or to the density functional method [8] for 
evaluating the van der Waals forces. This nonlinear 
formulation of an essentially linear problem simplifies 
the analysis of the effect of disorder in a similar way 
as it does in quantum problems [S]- We now write 
a,j(k, z) = ka^{k,z) and if the distributions of the 
ay'ik.z) — y are given by Pi(k,y) we find that, in three 
dimensions the average of the I dependent free energy is 
given by 

ksTA f f f 
(F) = — / dkk / dyi / dy 2 pi(k,y 1 )p 2 (k,y 2 ) 



4~ 



In 1 



(11) 



(yi + eo)(j/2 + eo) 

where the angled bracket on the l.h.s. indicates the dis- 
order average over the dielectric function within the slabs 
and we have assumed that the realizations of the disorder 
in the two slabs are independent. 

Let us first investigate the form of van der Waals in- 
teraction free energy in the limit of large separations be- 
tween the two slabs. The equation obeyed by a can be 
written as 

da(k, £) 



<C/kV 



(12) 



with £ = zk. When k is small e(£/&) varies very rapidly 
and so becomes decorrelated from the value of a. The 
Laplace transform for the probability density function of 
a is defined by p{k,s,() = J °° dyexp(-sy)p(k, y, C) = 
(exp(— sa(k, C))) and, from the equation of motion Eq. 



(121, obeys 



1 dp(k, s, C) = 
's d( 

e(C/k)exp(-sa(k,0) 



<C/k) 



exp(— sa{k, Q) 



Assuming that k is small and thus that a(k, Q and e((/k) 
are decorrelated, we can write 



1 dp(k, s, C) 
s d( 



■>',/,(/,. s.o- l,e)^(fc, S ,C). (14) 



As we are interested in the limit of thick slabs it suffices to 
know the equilibrium distribution of this equation which 
is given by lim^^ 00 p(fc, s, £) = cxp(— e*s) with 



(15) 



Inverting the Laplace transform then gives the equilib- 
rium distribution p(y, k) = 5{y — (*) at small k. When £ 
is large the integral in Eq. (Ill is dominated by the small 




k behavior and we may use the analysis presented above, 
to give the following asymptotic form for the interaction 
free energy 



(F)(e 



oo 



(16) 

with A* = (e* - eo)/(e| + £o) and where e* are defined 
via Eq. ( 15 1 . The subscript i on the angled brackets 



signifies that we are averaging the dielectric function in 
the slab i. The term H* defines an effective disorder- 
dependent Hamaker coefficient. This therefore justifies 
physical arguments replacing the random layered mate- 
rial by an effective anisotropic medium where the dielec- 



tric tensor is has the form e 



(e) _ 



and Si! 



M _ 



all other terms being zero by symmetry. The term 
£|| is the effective dielectric function in the z direction 

e[, e = l/(l/e), and the perpendicular components are 



II 

given by e 



(e) 



(e). The expressions for ejf' and e 



ie) 



follow simply from the fact that in the perpendicular 
direction the dielectric function is obtained by analogy 
to capacitors in series and in the parallel direction by 
analogy to capacitors in parallel arrangement [10 . The 
effective value, e* , for dielectric constant of this system 
coincides with that of Eq. (15) above [7]. This result 
shows that for large separations (where i is much larger 
than the correlation length of the dielectric disorder) the 
thermal Casimir interaction free energy is self averaging 
and agrees with that given by physical reasoning. 

One would imagine that as the distance between the 
slabs is reduced, the result will be increasingly dominated 
by the slab composition at the two opposite faces 
Indeed in the small £ limit Eq. (11) is dominated by 
the large k behavior. The asymptotic behavior can be 
extracted if one assumes the ansatz 



a(z,k) = 



71=0 



g n (z) 

k n 



(17) 



(13) Substituting this into Eq. (12 1 gives the following chain 
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of equations for a n (z) 
1^1 da n {z) _ 1 -A a n (z)a m (z) 

h 1^ h n dz ~ £ ( z ) 1^ k m+n ' \ 

n,m—0 



k ^ k r ' 

n=0 



From here it is easy to see that to order O(l) the leading 
asymptotic result of Eq. ( 23 1 is given by 



a (z) = e(z) 



(19) 



The equation for the corrections (n > 1) to this asymp- 
totic limit is 



da n _i(z) 
dz 



1 " 



dz) ^- 

v ' m=0 



and the next two terms from this expansion yield 
1 de(z) 



ai(z) 
a 2 (z) 



2 dz 
1 d 2 e(z) 
4 dz 2 8e(z) \ dz 



1 



de(z) 



(20) 

(21) 
(22) 



It is straightforward to realize that these terms generate 
0(1/ 1) corrections to the asymptotic result which are 
subdominant when £ is large. Thus to the leading order 

a(z, k) ~ oto(z) — e(z) (23) 

and from here it follows straightforwardly that 



lim Pi(y,k) = pi(e) 



(24) 



where pi is the probability density function of e(z) in 
medium i. This result is easily understood from the phys- 
ical discussion above. The average of the thermal Casimir 



interaction free energy Eq. (Ill in the small separation 
limit is then given by 



k B TA 



idu 



I6ir£ 2 

In (1 - AiA 2 e~ u 



Pi{ti)P2{£2)deide 2 

(25) 



with A,; = (cj — eo)/(£j + e o)- The forms of the thermal 
Casimir interaction free energy are thus given by Eqs. 



(25 1 and (16 1 in the small and large interslab separation 



limits respectively. 

In the limit of large separation between the slabs 
we have obtained the limiting behavior of the thermal 
Casimir effect and shown that the free energy is given 
by self-averaging and that the distributions of a(k, z) are 
strongly peaked. It can be shown [7] that the attraction 
at large separation between two (statistically identical) 
homogeneous media (with e = (e)) is stronger than that 
between the two fluctuating media if (l/e) -1 > eo. How- 
ever it is always weaker if (e) < eo- So, depending on the 
details of the distribution of the fluctuating dielectric re- 
sponse in the two slabs and the dielectric response of 



the medium in-between, the effective interaction at large 
inter-slab separations can be stronger or weaker than that 
for a uniform medium with a dielectric constant equal to 
the mean dielectric function of the fluctuating media. 

For small separations the interaction free energy is a 
random variable which has to be averaged over the prob- 
ability density function of the dielectric functions in the 
media composing the two interacting slabs. The interme- 
diate length scales can be analyzed via perturbation the- 
ory [7] , and there may also exist models of disorder that 
can be treated exactly. The nonlinear formulation of the 
problem presented here should be equally useful to treat 
the case of deterministically varying dielectric functions 
and could open up a useful computational framework for 
designing materials where the effective interaction can be 
tuned to induce attractive or repulsive forces depending 
on the separation, for practical applications |llj . 
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